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We study spherically symmetric static spacetimes generally filled with an anisotropic fluid in the 
nonrelativistic general covariant theory of gravity. In particular, we find that the vacuum solutions 
are not unique, and can be expressed in terms of the U(l) gauge field A. When solar system tests 
are considered, severe constraints on A are obtained, which seemingly pick up the Schwarzschild 
solution uniquely. In contrast to other versions of the Horava-Lifshitz theory, non-singular static 
stars made of a perfect fluid without heat flow can be constructed, due to the coupling of the fluid 
with the gauge field. These include the solutions with a constant pressure. We also study the general 
junction conditions across the surface of a star. In general, the conditions allow the existence of a 
thin matter shell on the surface. When applying these conditions to the perfect fluid solutions with 
the vacuum ones as describing their external spacetimes, we find explicitly the matching conditions 
in terms of the parameters appearing in the solutions. Such matching is possible even without the 
presence of a thin matter shell. 
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I. INTRODUCTION 

Recently, Horava proposed a theory of quantum grav- 
ity [l[, motivated by the Lifshitz theory in solid state 
physics @. Due to several remarkable features, the 
Horava-Lifshitz (HL) theory has attracted a great deal 
of attention (see for example, Q and references therein) . 

In the HL theory, Lorentz symmetry is broken in the 
ultraviolet (UV). The breaking manifests in the strong 
anisotropic scalings between space and time, 

X ->•&£, t^tlH. (1.1) 

In (3 + l)-dimensional spacetimes, the theory is power- 
counting renormalizable, provided that z > 3. At low 
energies, the theory is expected to flow to z — 1, whereby 
the Lorentz invariance is "accidentally restored." Such an 
anisotropy between time and space can be easily realized, 
when writing the metric in the Arnowitt-Deser-Misner 
(ADM) form 

ds 2 = -N 2 c 2 dt 2 + gij (dx l + N l dt) {dx J + N j dt) , 

(i, j = 1,2,3). (1.2) 

Under the rescaling (jl.ip with z = 3 (a condition we shall 
assume in this paper), N, N l and g^ scale as, 

N N, iV* — > t~ 2 N\ 9ij ->■ 9ij . (1.3) 

The gauge symmetry of the system are the foliation- 
preserving diffeomorphisms Diff(M, J 7 ), 

t = t-f(t), ^ = ^-C(t,x), (1.4) 
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for which N, N z and transform as 

5 gij = ViCj + VjO + ./.'/•.; • 

6N t = N k Vi£ k + t k VkNi+g ik <; k + Nif + Nif, 

SN = ( k V k N + Nf + Nf, (1.5) 

where / = df/dt, V; denotes the covariant derivative 
with respect to the 3-metric g^, and iVj = gi k N k , etc. 
From these expressions one can see that the lapse func- 
tion N and the shift vector iVj play the role of gauge fields 
of the Diff(M, F) symmetry. Therefore, it is natural to 
assume that iV and Ni inherit the same dependence on 
spacetime as the corresponding generators, 

N = N(t), Ni = Ni(t,x), (1.6) 

which is clearly preserved by the Diff (M, F) , and usually 
referred to as the projectability condition. 

Due to the restricted diffeomorphisms (|1.4|) . one more 
degree of freedom appears in the gravitational sector - the 
spin-0 graviton. This is potentially dangerous, and needs 
to decouple in the infrared (IR) , in order to be consistent 
with observations. Whether it is the case or not is still 
an open question [f|. In particular, the spin-0 mode is 
not stable in the Minkowski background, in the original 
version of the HL theory and in the Sotiriou, Visser 
and Weinfurtner (SVW) generalization @, 0j- Although 
in the SVW setup it is stable in the de Sitter background 
[|| . In addition, non-perturbative analysis showed that it 
indeed decouples in the vacuum spherical static [j| and 
cosmological [8| spacetimes. 

To overcome the problem, various models have been 
proposed Q. In particular, Horava and Melby- 
Thompson (HMT) [10| recently proposed a version in 
which the spin-0 graviton is completely eliminated by 
introducing a Newtonian pre-potential ip and a local 
U(l) gauge field A, so that the foliation-preserving- 
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diffeomorphisms, Diff(M, 7), are extended to 

17(1) k Diff(M, 7). (1.7) 

Effectively, the spatial diffeomorphism symmetries of 
general relativity are kept intact, but the time 
reparametrization symmetry is contracted to a local 
gauge symmetry [11] . The restoration of general covari- 
ance, characterized by Eq. (jl.7p . nicely maintains the spe- 
cial status of time, so that the anisotropic scaling (II. ip 
with z > 1 can still be realized. Under the Diff(A7, 7), 
the fields A and if transform as, 

5 A = (%A + fA + fA, 

Sip = ftp + CBiV, (1-8) 

while under the local [7(1), they, together with g^, trans- 
form as 

S a A = a — Af'Via, J Q <^ = —a, 
S a Ni = NVia, 5 agij = = S a N, (1.9) 

where a is the generator of the local [7(1) gauge symme- 
try. A remarkable by-production of this "non-relativistic 
general covariant" setup is that it forces the coupling 
constant A, introduced originally to characterize the de- 
viation of the kinetic part of the action from GR [l| , to 
take exactly its relativistic value A = 1. A different view 
can be found in [l2j | . 

In this paper, we investigate systematically the spher- 
ically symmetric spacetimes in the HMT setup. In par- 
ticular, after briefly reviewing the theory in Sec. II, we 
develop the general formulas of such spacetimes gener- 
ally filled with an anisotropic fluid with heat flow in Sec. 
III. Then, in Sec. IV we study the vacuum solutions, 
and express them all in terms of the gauge field A. Al- 
though the solutions are not unique, when we apply them 
to solar system tests in Sec. V, we find that these tests 
seemingly pick up the Schwarzschild solution generically. 
It should be noted that solar system tests were studied 
by several authors in the framework of the HL theory 
mainly without the projectability condition [l3| . In Sec. 
VI, we study perfect fluid solutions without heat flow, 
and also express them all in terms of the gauge field A. 
By properly choosing A, non-singular stars can be con- 
structed, due to the coupling of the gauge field with the 
fluid. This is in contrast to all the previous versions of 
the HL theory, in which it was shown that non-singular 
sta tic p erfect fluid soltuons without heat flow do not ex- 
ist [14| . although the ones with heat fluid do [l6j]. Then, 
we restrict ourselves to the cases where the pressure is a 
constant, which is quite similar to the Schwarzshcild fluid 
solution. We show explicitly that these solutions are free 
of spacetime singularities at the center of the star. In Sec. 
VII we consider the junction conditions across the surface 
of a compact object with the minimal requirement that 
the matching is mathematically meaningful. In particu- 
lar, this allows the existence of a thin matter shell on the 
surface of the star, where the pressures of the thin shell 



can have high-order derivatives of the Dirac delta func- 
tion, applying these conditions to the solutions found in 
Sees. V and VI, we obtain the matching conditions with 
or without a thin matter shell. Finally, in Sec. VIII we 
present our main results and concluding remarks. 

It should be noted that spherically symetric static 
spacetimes in other versions of the HL theory have been 
extensively studied, and it is difficult to give a complete 
list here. Instead, we simply refer readers to references 
given in Q , and to the more recent review articles [!, 0] . 

II. NONRELATIVISITC GENERAL 
COVARIANT HL THEORY 

In this section, we shall give a very brief introduction 
to the non-relativistic general covariant theory of gravity. 
For details, we refer readers to [HJ. We shall closely 
follow [15J, so that the notations and conversations will 
be used directly without further explanations. In the 
following, [l5[ will be referred to as Paper I. 

The total action is given by, 

S = Q 2 J dtd 3 xN^/g(£ K -£ v +£ tp + £ A 

+^£m) , (2.1) 

where g = det gij , and 

£ K = KijK 13 — K 2 , 

C v = ,;t/"(2A,, • V,V,,::). 

C A = ^( 2 A s -i?). (2.2) 

Here the coupling constant A g , acting like a 3- 
dimensional cosmological constant, has the dimension of 
(length) ~ 2 . The Ricci and Riemann terms all refer to the 
three-metric g^. Kij is the extrinsic curvature, and Gij is 
the 3-dimensional "generalized" Einstein tensor, defined 
by 

Gij = />',, ^!h,l> I \, !<!,,■ (2.3) 

£m is the matter Lagrangian density and £y an arbitrary 
Diff(S)-invariant local scalar functional built out of the 
spatial metric, its Riemann tensor and spatial covariant 
derivatives, without the use of time derivatives. In [6], by 
assuming that the highest order derivatives are six and 
that the theory respects the parity, SVW constructed the 
most general form of £y, given by 

Cv = ( 2 9o + giR+^(92R 2 + g3RijR ij ) 
+^ (g- if + //:,/.' RijR ij + //,//;//; //:) 

+ ^ [g 7 RV 2 R + g 8 (V,R jk ) (V l W k )] , (2.4) 



3 



where the coupling constants g s (s = 0, 1, 2, . . . 8) are all 
dimensionless. The relativistic limit in the IR requires 
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■i, c 2 



1 



16ttG 



(2.5) 



In Paper I, this possibility was left open. To compare 
with the results obtained in (T^ |. which will be referred 
to as Paper II, we shall restrict ourselves to these values. 

Variation of the total action (|2.ip with respect to the 
lapse function N(t) yields the Hamiltonian constraint, 

J &xy/g {L K + £v - v& s ViVj<p) = 8nG J d?x^gj\ 

(2.6) 

where 



n _J(NC M ) 



J = 2- 



SN 



(2.7) 



Variation of the action with respect to the shift N* 
yields the super-momentum constraint, 



V ; (-'' ,:V'") = 8nGJ\ 



(2.8) 



where the super-momentum Tt lJ and matter current J 1 
are defined as 



SC 



K 



l\" ■ Kg' 1 . 



J' 



-N 



(2.9) 



Similarly, variations of the action with respect to <p and 
A yield, 



G ij (My + V, V^) = 8ttGJ l 
R-2A g = 8ttGJ a , 



where 



SC 



M 



8 ip 



Ja 



S(NC 



M) 



5A 



(2.10) 
(2.11) 

(2.12) 



On the other hand, variation with respect to gij leads to 
the dynamical equations, 



1 



N^g 



[y/g {it 13 - yQ ij )] = -2 (K 2 ) 13 + 2KK i3 



N 
1 



:V fe 



N k ir lj - 2tt* (< JV j) 



+ F» + F l J + Fjl + SitGtV , (2.13) 
where (K 2 ) ij = K a Kf, f m = (/ y + fa) /2, and 

l r 



r A 



N 



AR ij - V¥-3 ij V 2 4 



<p 



n=l 



F 13 = 



i s (-wvl = £ gaC . {Fa y 



(2.14) 



with n s = (2,0,-2,-2,-4,-4,-4,-4,-4). The 
geometric 3-tensors (Fg)^ and Ffe n ) are given by 
Eqs.(2.21)-(2.23) in Paper I, which, for the sake of the 
readers' convenience, are reproduced in Eqs. (|A.l[) and 
(IA.2I) of this paper. The stress 3-tensor r %3 is defined as 



2 S(^gC M ) 



(2.15) 



The matter quantities (J*, J\ J v , Ja, t %3 ) satisfy 
the conservation laws, 



J d 3 xyfg 



(V9Ja) 



0, (2.16) 



Ni 
' N 



V k J k + J^iip - ^ViA = 0. (2.17) 



III. SPHERICALLY SYMMETRIC STATIC 
SPACETIMES 

Spherically symmetric static spacetimes in the frame- 
work of the SVW setup are studied systematically in Pa- 
per II. In this section, we shall closely follow the develop- 
ment presented there. In particular, the metric for static 
spherically symmetric spacetimes that preserve the form 
of Eq. (|1.2[) with the projectability condition can be cast 
in the form, 



ds 2 



c 2 dt 2 + e 2lJ (dr + e^dt) + r 2 d 2 tt, (3.1) 



in the spherical coordinates x 1 — (r,8,4>), where d 2 il 



d9 2 + sin 2 9d6 2 , and 



V = n{r), v = y{r), JV* = {e"-', 0, 0} 



(3.2) 



The corresponding timclikc Killing vector is £ — dt- For 
the above metric, one finds 



2v' 



-2v 



T-e 



2i/\ 



a, 



c 



K 



4 



r — 



tpe 



if" - i/<p' 



( V 2 i) + 1] (» 

+2e"+V) - l[v' - A 9 re 2 ") (V + 2e^ +lJ ^j j, 
^ = ^[e-^(l-2rV) + (A 9 r 2 -l)' 
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0) 



(3.3) 



s=0 



where a prime denotes the ordinary derivative with re- 
spect to its indicated argument, = <5| + sin 2 OSfsj, 

and jC^'s are given by Eq.(Al) in Paper II. Then, the 
Hamiltonian constraint (I2.6[) reads, 

/" ^ ic + £v-£W-87rGJ*)e !/ r 2 dr = 0, (3.4) 



where 



(A g r 2 - 1) + 1 (p" - iV) 



Fij's for the metric (13. 1|) are given by Eq.(A2) in Paper 
II, and Iff s -, are given by Eg. (IB. II) in the present paper. 
As in Paper II, here we define a fluid with p r — pg as a 
perfect fluid, which in general conducts heat flow along 
the radial direction (l7j . 

Since the spacetime is static, one can see that now 
the energy conservation law (|2.16p is satisfied identically, 
while the momentum conservation (|2.17p yields, 

2 1 

v(i'-(v'-p' r )--(v-p r +pg)+J tp (p'--J A A' = 0. (3.12) 

To relate the quantities J*, J % and ry to the ones often 
used in general relativity, following Paper II, one can first 
introduce the unit normal vector rt M to the hypersurfaces 
t = Constant, and then the spacelike unit vectors x M , 9^ 
and <j>u, defined as [l8| 



3 



(3.13) 



(3.5) 



while the momentum constraint (|2.8|) yields, 

2rv' + e-^ +v) [e 2l/ (A g r 2 - l) + 1 if 

= -8nGr 2 e 2 ^-^v, (3.6) 

where J 1 = e^^ +v ' (y, 0, 0). It can be also shown that 
Eqs. (|2~TU|) and (|2~TTj) now read, 



e 2v (A g r 2 - 1) + 1 



ip — v ip + /i 



-2{v' - A a re 2vXl J 1 "" +v ^ - °~^~ 2 ^ v 



ip' + e^+n = %irGr 2 e^J Vl (3.7) 



2rv' 



,2w 



(A s r 2 - 1) + 1 



^Gr 2 e 2v J A 



(3.8) 



The dynamical equations (|2.13[) . on the other hand, yield, 

2(/x' + i/') + " + ^ e2( "" M) + 
r 2 

= -re" 2 " (i^ rr + + + 87rGe 2 V) , (3.9) 
//' + (2// - i/) (V + ±) + l -e 2 ^^ (C v + C A ) 



Fee + F£ g + F d A d + 8nGr 2 P e , (3.10) 



where 



= e 2 y r 5l8 r 3 + r 2 p e n 



ij • 



= -(A' + Av') + e~ 2v [r 2 (A" - v' A') 

+ r(A' + Av') - A(l- e 2,y )]n %3 , (3.11) 



In terms of these four unit vectors, the energy-momentum 
tensor for an anisotropic fluid with heat flow can be writ- 
ten as 

+PrXnXv + Pe (0»0v + <t>n<t>v) , (3-14) 

where pn, q, p r and pg denote, respectively, the en- 
ergy density, heat flow along radial direction, radial, and 
tangential pressures, measured by the observer with the 
four- velocity n^. Then, one can see that such a decompo- 
sition is consistent with the quantities J* and J 1 , defined 
by 



PH 



-2f 



(3.15) 



It should be noted that the definitions of the energy den- 
sity ph, the radial pressure p r and the heat flow q are 
different from the ones defined in a comoving frame in 
general relativity. For detail, we refer readers to Ap- 
pendix B of Paper II. 

Finally, we note that in writing the above equations, we 
leave the choice of the U(l) gauge open. From Eq. (|1.9[) 
one can see that it can be used to set one (and only one) 
of the three functions A, ip and N r to zero. To compare 
our results with the one obtained in [10|, in the rest of 
this paper (except the first part of Sec. VII), without 
loss of the generality, we shall choose the gauge, 



p = 0. 

Then, we find that 

A= = 0, F{l, n) =0, (n = l,2,3). 



(3.16) 



(3.17) 
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IV. VACUUM SOLUTIONS 



SOLAR SYSTEM TESTS 



In the vacuum case, we have J* = v = p r = pg = J a = 
J v = 0. With the gauge (|3.16[) . from the momentum 
constraint (|3.6[) we immediately obtain v = Constant, 
while Eq. (|3.8I) further requires 



V = 0, kg = 0. (4.1) 

This is different from the solutions presented in [Icj . 
where ^ ^ 0, /i = — oo. Inserting the above into Eq. (|3.7p . 
it can be shown that it is satisfied identically. Since 
v = 0, from the expressions of (-F s )y given by Eq.(A2) 
in Paper II, we find that = for s ^ 0, and 

(-Fb)y = -9ij/2, so that 

F« = -A 9ij , (4.2) 

where A = C 2 5o/ 2 - Substituting Eqs. (|3~TTj) and (|3~T7)) - 
(pi~2jl into Eqs.dHJ]) and (f3TT0|) . we find that 



(2r/i' + l)( 



,2/- 



Ar 2 



2//U' + - 



- 2rA', 

-2ju 



Ar- (rvi')' 



(4.3) 
(4.4) 



It can be shown that Eq. (|4.4l) is not independent, and 
can be obtained from Eq. (|4.3[) . Therefore, the solutions 
are not uniquely determined, since now we have only one 
equation, (|4.3p . for two unknowns, fi and A. In particu- 
lar, for any given A, from Eq. (|4.3[) we find that 



2m 1 , 
h -Ar 



2A{r) + ^- J A(r')dr' 



(4.5) 



On the other hand, we also have 
4A' 

C K = 2A, £ y = 2A. 

r 

Inserting it into the Hamiltonian constraint (I3.4p , we find 
that 



(4.6) 



rA'(r)dr = 0. 



(4.7) 



Therefore, for any given function A, subjected to the 
above constraint, the solutions given by Eqs. (|4.1[ ) and 
(|4.5[) represent the vacuum solutions of the HL theory. 
Thus, in contrast to general relativity, the vacuum solu- 
tions in the HMT setup are not unique. 

When A is a constant (without loss of generality, we 
can set A = 0), from the above we find that 



1, /2m 
/i= 2 hl 



+ ^Ar 2 ) , (A = 0), (4.8) 



which is exactly the Schwarzschild (anti-) de Sitter solu- 
tion, written in the Gullstrand-Painleve coordinates [lj| . 
It is interesting to note that when m = we must assume 
that A > 0, in order to have fi real. That is, the anti-de 
Sitter solution cannot be written in the static ADM form 

ZD. 



The solar system tests are usually written in terms 
of the Eddington parameters, by following the so-called 
"parameterized post-Newtonian" (PPN) approach, intro- 
duced initially by Eddington [20]. The gravitational field, 
produced by a point-like and motion-less particle with 
mass M, is often described by the form of metric (21.1, 



ds 2 = 



^c 2 dt 2 



2* 



dr 2 



'■d 2 n, 



(5-1) 



where $ and $ are functions of the dimensionless quan- 
tity x = GM/(rc 2 ) only. For the solar system, we have 
GMq/c 2 ~ 1.5km, so that in most cases we have x <C 1. 
Expanding ^ and $ in terms of x, we have [2l| 



„2* 



2$ 



1 



l + 2 7 



\ c z r J 



„,„ v fGMY 



(5.2) 



where /3 and 7 are the Eddington parameters. General 
relativity predicts /3 — 1 = 7 strictly, while the current 
radar ranging of the Cassini pro be 1 221) . and the proces- 
sion of lunar laser ranging data [23j yield, respectively, 
the bounds [H, 



7 - 1 = (2.1 ±2.3) x 10" 
j3- 1 = (1.2 ± 1.1) x 10" 



(5.3) 



which are consistent with the predictions of general rel- 
ativity. 

To apply the solar system tests to the HL theory, we 
need first to transfer the above bounds to the metric coef- 
ficients /j, and v. In Appendix B of Paper II, the relations 
between ($, ^) and (/x, v) have been worked out explic- 
itly, and are given by 



1 



/' : 



In 



(l-e-) 



c 2 (l 



3> + 



(5.4) 



or inversely, 



$ = v - - In ( 1 - -U 2m 
2 V c l 

^ = il n fi_i e 2M 



(5.5) 



Inserting Eq. (|5.2[) into Eq. (|5.4p . we find that 



11 = 2 ln< ^ 2c2 



( c 2 r ) t) ^ c 2 r ^ 

, ^ (GM\ 
(7-1) ^ +•••• 



(5.6) 



Comparing Eq. (|5.6p with Eq. (|4.5l) for A = 0, we find that 
in order to be consistent with solar system tests, we must 
assume that 



A(r) = O 



\ c z r J 



(5.7) 
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Together with the Hamiltonian constraint (|4.7I) . we find 
that this is impossible unless A = 0. Therefore, although 
the vacuum solution in the HMT setup is not unique, the 
solar system tests seemingly require that it must be the 
Schwarzschild vacuum solution. 

It should be noted that by choosing A(r) in very par- 
ticular forms, the condition A = could be relaxed [25| . 
But, such chosen A is not analytic (in terms of the dimen- 
sionless quantity and it is not clear how to expand 
it in the form of (|5.6p . Thus, in this paper we simply 
discard those possibilities. 



Substituting it into Eq. (|6.6[) . and then taking a deriva- 
tive of it, we find that the resulting equation is exactly 
given by Eq. ([6l8]) . Thus, both Eqs. (g5§ and flTTf are 
not independent, and can all be derived from Eqs. ([6.6[ ) 
and (I6.10p . Then, in the present case there are five in- 
dependent equations, the Hamiltonian constraint (13.41) . 
and Eqs. (|6.4[) . (I6.5[) . (|6.6p and (|6.10l) . However, we have 
six unknowns, A, p, p, J*, J v , J a- Therefore, the prob- 
lem now is not uniquely determined. As in the vacuum 
case, we can express all these quantities in terms of the 
gauge field A. In particular, substituting Eq. (|6.10[) into 
Eq. ([6.6[) and then integrating it, we obtain 



VI. PERFECT FLUID SOLUTIONS 

In this section, let us consider perfect fluid without 
heat flow, that is, 



Vr = Pe =P, v = 0. 



(6.1) 



Then, together with the gauge choice (|3 . 1 6|) . from 
Eq. (|3.6[) we find that v = Constant. However, to be 
matched with the vacuum solutions outside of the fluid, 
as shown in Sec. IV, we must set this constant to zero, 



v = 0, 



(6.2) 



from which we immediately find that Rij = 0, is still 
given by Eq. (|4.2[) , and 



C A = 2A g A, 



F, 



2A' 



(6.3) 



Inserting the above into Eqs. p.7p - (|3.10[) , we find that 
A 



(6.4) 
(6.5) 



(r/)' + 2rA' + A g r 2 A - Ar 2 = -8irGr 2 p, (6.6) 
\rf + f+ (rA')' + A g rA - Ar = SirGrp, (6.7) 

where / = e 2/i . From the last two equations, we find that 

r 2 /"-2/ = -2r 3 

On the other hand, the conservation law of momentum 
(|3.12p now reduces to 



P 



8vrG 



A' = 0, 



which has the solution, 



P = Po 



8ttG 



A. 



(6.9) 



(6.10) 



where po is an integration constant. 



u = — In 
p 2 



— + -(A-87TGp Q )r 2 -2A+- f A(r')dr 
r 3 V / r J 



3.11) 



Then, we find that 

C v = 2(A-87rGp ), C v = 2A. (6.12) 

r 

Inserting the above into Eq. (l3.4p . we find that it can be 
cast in the form, 



p{r)dr = 0, 



where 



2ttG 



(6.13) 



3.14) 



From the above one can see that once A is given, one 
can immediately obtain all the rest. By properly choos- 
ing it (and p(r)), it is not difficult to see that one can 
construct non-singular solutions representing stars made 
of a perfect fluid without heat flow. To see this explicitly, 
let us consider the following two particular cases. 



A =0 



When A g — 0, we have 



J v = Ja = 0, p = Po, 



(6.15) 



while p and J* are still given by Eqs. (|6.1ip and (|6 . 14[) . 
respectively. To have a physically acceptable model, we 
require that the fluid be non-singular in the center. Since 
Rij = 0, one can see that any quantity built from the 
Riemann and Ricci tensor vanishes in the present case. 
Then, possible singularities can only come from the ki- 
netic part, Kij, where the very first quantity is 



K 



<r'K,j = - (r/x' + 2) 
r 

e —— f — + (A - 8irGpo) r 2 - rA' - 3A 
r \ r 



+- I A{r')dr' 



(6.16) 
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Assuming that near the center A is dominated by the 
term r a , we find that K is non-singular only when 



m = 0, a > 2. 



3.17) 



For such a function A, Eq. ([6.14[) show that J* is non- 
simgular, as long as p(r) ~ 0(r 2 ). 



B. A = A 



When A is a constant, from Eq. (|6.10|) we can see that 
the pressure p is also a constant. Then, the integration 
of Eq.HU) yields, 




A - 8ttGpq r 



(6.18) 



Inserting it into Eq. (|6.8[) we find that it is satisfied iden- 
tically, while the Hamiltonian constraint (|3.4I) can also 
be cast in the form of Eq. (|6.13[) , but now with 



J* = 



1 (lOKGpo-eQ 



8nG V 



(6.19) 



Thus, the solutions of Eqs. (|6.2[) and (|6.18[) represent a 
perfect fluid with a constant pressure, p — p(Aq). In this 
case, it can be shown that K is free of any spacetime 
singularity at the center only when m = 0. 

It should be noted that in [T3| it was shown that non- 
singular static solutions of perfect fluid without heat flow 
do not exist. Since their conclusions only come from the 
conservation law of momentum, one might expect that 
this is also true in the current setup. However, from 
Eq. (l3.12l) we can see that in the present case the conser- 
vation law contains two extra terms, J v and J a- Only 
when both of them vanish, can one obtain the above con- 
clusions. Since in general one can only choose one of them 
to be zero by using the gauge freedom, it is expected that 
non-singular static stars can be constrcuted by properly 
choosing the gauge field A. 

It should be also noted that the arguments presented in 
14] do not apply to the case where the pressure is a con- 
stant. Therefore, when p' = non-singular stars without 
heat flow can be also constrcuted in other versions of the 
HL theory, although when p 1 ^ this is possible only in 
the HMT setup. In addition, the definitions of the quan- 
tities p, p and v (= qe^) used in this paper are different 
from the ones used usually in general relativity. For de- 
tail, we refer readers to Appendix B of [l6|, specially to 
Eq.(B16). 



VII. JUNCTION CONDITIONS 

To consider the junction conditions across the hyper- 
surface of a compact object, let us first divide the whole 
spacetime into three regions, and S, where V~ (V + ) 



denotes the internal (external) region of the star, and 
E is the surface of the star. As shown in Paper II, 
once the metric is cast in the form (13.11) . the coordi- 
nates t and r are all uniquely defined, so that the co- 
ordinates defined in V + and V~ must be the same, 
= {a;"^} = (t,r,9,4>). Since the quadratic terms 
of the highest derivatives of the metric coefficients p, and 
v are only terms of the forms, v" 2 , v"v"' and // 2 , the 
minimal requirements for these two functions are that 
u(r) and p(r) are respectively at least C 1 and C° across 
the surface E, and that they are at least C 4 and C 1 else- 
where. For detail, we refer readers to (l6| . Similarly, the 
quadratic terms of the Newtonian pre-potential are only 
involved with the forms, (pip", p 2 r , and (pip'. Therefore, 
the minimal requirment for ip is to be at least C° across 
the surface S. On the other hand, the gauge field A and 
its derivatives all appear linearly. Thus, mathematically 
it can be even not continuous across E. However, in this 
paper we shall require that A be at least C° too across 
E. Elsewhere, A and ip are at least C 1 . Then, we can 
write A and ip in the form, 

E{r) = E+(r)H (r - r ) + E~ (r) [1 - H (r - r )} , 

where E = [A, tp), ro is the radius of the star, and H (x) 
denotes the Heavside function, defined as 



H(x) 



(I, x>0, 
10, x<0. 



Since A and <p are continuous (C°) across 
must have 



(7.2) 



r , we 



]imit r ^ r +E+(r) = limits- E~(r). (7.3) 
Then, we find that 

E'(r) = E D r (r), 

E"(r) - E D rr (r) + [E']-S(r-r ), (7.4) 

where 

[E'Y = hmit r ^,+ £+(r) -limits- E~(r), 
E D {r) = E + H(r-r ) + E-\l-H(r-r j\.(7.5) 

Combining the above with Eq.(6.6) of Paper II, we find 
that 



l-'K — t-'K-i — l~Ai 



(7.6) 



where 



C 



Irn 



8g 7 e 



-6v 



C 



I in 



pe 



C 4 r 3 



2rv' 



k g r 2 e 2v + (1 - e 2u ) [<p'}. (7.7) 
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Setting 



J = J D + J Im 6(r - r ), 



(7- 



where J = { J*, u, J v , J a}, and J has support only 
on E, we find that the Hamiltonian constraint (|3.4j) can 
be written as 



C K +C v - tfP - 8ttGJ* ) e v r 2 dr 



471" V * 



V h 



(7.9) 



where 

cD 



(rro—e roo 
/ I(r)dr + / I(r)dr 
JO Jr +e 

(7.10) 

It shold be noted that in writing Eq. (j7.9[) , we had used 
the conversion, 



^d 3 xf(r)5(r - r ) = /(r ). 



(7.11) 



The momentum constraint (|3.6p will take the same 
form in Regions V^, while on the surface E it yields 

v Im = 0. (7.12) 

That is, the surface does not support impulsive heat flow 
in the radial direction. Similarly, Eqs. (|3.7p and (|3.8|) take 
the same forms in Regions . While on E they reduce, 
respectively, to 



[if/] = 8TrGr 2 e^J^ m , (7.13) 
(7.14) 



A g r 2 e 2v + (1 - e 21 
J A m = 0, (r = r ). 

On the other hand, the dynamical equations Q3.9[) and 
(I3.10[) take the same forms in Regions V^, and on the 
surface E they yield, 



-2v 



A g r 2 e 2 » + (l~e 2 ») [cp' 



+ 2r 2 F% Im y{r-r a ) 
-2r 2 (f/™ + 8nGe 2u p I r m ) , (7.15) 



| [fJ-T + \e 2 ^CJ™ 

{F^ m + F e A / m )\s(r-r ) 



(F^ + SnGr 2 ^), (7.16) 



where FL™ are given by Eq.(A5) in Paper II, Ff-' Im are 



given by Eq. (IC.2[ ) in Appendix C of this paper, and 

L = L D + L Im : (7.17) 



with L = (p r , pe). From Eq.(A5) in Paper II we can see 
that L Im in general takes the form, 

L lm = LW Im S(r-ro)+L^ Im 8'(r-ro)+L^ Im 8''(r-r ). 

The above represents the general junction conditions of a 
spherical compact object made of a fluid with heat flow, 
in which a thin matter shell appears on E. 

In the following we shall consider the matching of the 
perfect fluid solutions found in Sec. VI to the vacuum 
ones found in Sec. V. Since 



v ± = 0, ip ± 



0. 



(7.19) 



0, ££ = 2A±, £\T = 0, 



we immediately obtain 

4 - " r 

F£ = -A ±5 ±, F^ = 0, 
C% = C 1 ™ = 0, .. = (F*™) .. = 0. (7.20) 

Then, from Eqs.JEH, fEUl, (|7T4)l and (j7TT5l) we find 



,7m 



j 



Ira 



J 



Im 
A 



■Vr m 



0. 



(7.21) 



That is, the radial pressure of the thin shell must vanish. 
This is similar to what happened in the relativistic case 
17]. _ 

In the external region, V + , the spacetime is vacuum, 
and the general solutions are given by Eq. (14.5[) . 



1 



hi 



— + -A+r 2 -2A+ r 
r 3 



for which we have 



4A+ 



- 2A. 



£+=0. 



A + {r')dr' 

(7.22) 

(7.23) 



In the internal region, V~, two classes of solutions of 
perfect fluid without heat flow are found, given, respec- 
tively, by Eqs. (|6.11l) and (|6.18l) . which can be written 



as, 



A_ -8nGp )r z -2A~{r) 



+ - / A-{r')dr' \, 



for A g — 0, and 



p = i In < ^ A_ - A- A - 8irGp 



(7.24) 



(7.25) 



for A — Aq, where p = A „A 2 /2 +po- Then, we find 



that 



c - = I 2(87rGpo - A_) + —f-, A" = 0, 

K \ 2(8ttG P - A_ + A-^o), A-=A , 
C A = 2A~A-(r). (7.26) 

To further study the junction conditions, let us consider 
the two cases A~ = and A~ — Aq separately. 
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A. Ag = 



In this case, the continuity conditions of /x and A across 
£ read, 



1 r r o 
m + -AArg - 

A+(r ) = A-(r ), 



AA(r)dr 



AitG 



-Por'o, 



(7.27) 



where AA = A + - A_ and AA = A + - A~ . Then, the 
Hamiltonian constraint (I7.9[) becomes, 

p(r)dr+ rA+(r)dr = -GJ tJm , (7.28) 

Jro ^ 

while the dynamical equation (|7.16[) reduces to, 

AA = -8ttG{po + 2^ 0)/m ) , (7.29) 

where y\ m = pf )Im 5{r - r ) [cf. Ea. fTTgl) ]. 

When the matter thin shell does not exist, we must 
set J*< /m = p ( g )Im = 0, and Eos. (17371) - (17391) become 
the matching conditions for the constants A±, m, pg 
and the functions A^(r) and p(r). 

B. A = A Q 

In this case, it can be shown that the continuity con- 
ditions for fi and A become, 

1 f r ° 
m + - AArg + / A+(r)dr = r A Q 

- ^(A-A + 8TTGp)rf ) , 
A+(r )=A Q , (7.30) 
while the Hamiltonian constraint (|7.9j) reduces to, 

i roo 

p{r)dr +4 / rA+(r)rfr = 2GJ tIm . (7.31) 
The dynamical equation (|7.16|) . on the other hand, yields, 



AA = -A- A 



(7.32) 



In all the above cases, one can see that the matching is 
possible even without a thin matter shell on the surface 
of the star, J*> /m = = p^ Im ', by properly choosing the 
free parameters. 



VIII. CONCLUSIONS 

In this paper, we have systematically studied spher- 
ically symmetric static spacetimes generally filled with 
an anisotropic fluid with heat flow along the radial direc- 
tion. When the spacetimes are vacuum, we have found 



solutions, given explicitly by Eqs. (|4.1[) and (14.5[) . from 
which one can see that the solution is not unique, be- 
cause the gauge field A is undetermined. When A = 0, 
the solutions reduce to the Schwarzschild (anti-) de Sit- 
ter solution. We have also studied the solar system tests, 
and found the constraint on the choice of A. 

It should be noted that we have adopted a different 
point of view of the gauge field A in the IR limit than 
that adopted in [T3|. In this paper we have considered 
it as independent from the 4D metric g^ u , although it 
interacts with them through the field equations. This 
is quite similar to the Brans-Dicke (BD) scalar field in 
the BD theory, where the scalar field represents a de- 
gree of freedom of gravity, is independent of the metric, 
and its effects to the spacetime are only through the field 
equations 26]. On the contrary, in [l0| the authors con- 
sidered the gauge field A as a part of the lapse function, 
g tt ~ — {N — A) 2 in the IR limit. 

We have also investigated anisotropic fluids with 
heat flow, and found perfect fluid solutions, given by 
Eq. ([6. lip . By properly choosing the gauge field A, the so- 
lutions can be free of spacetime singularities at the center. 
This is in contrast to other versions oftheHL theory 
due to the coupling of the fluid with the gauge field. We 
then have considered two particular cases, in which the 
pressure is a constant, quite similar to the Schwarzschild 
perfect fluid solution. In all these cases, the spacetimes 
are free of singularities at the center. 

For a compact object, the spacetime outside of it is 
vacuum, matching conditions are needed across the sur- 
face of the star. With the minimal requirement that the 
junctions be mathematically meaningful, we have worked 
out the general matching conditions, given by Eqs. (17.9[) 
and (|7.12[) - (|7.16j) . in which a thin matter shell in gen- 
eral appears on the surface of the star. Applying them 
to the perfect fluids, where the spacetime outside is de- 
scribed by the vacuum solutions (|4.5p , we have found the 
matching conditions in terms of the free parameters of 
the solutions. When the thin shell is removed, these con- 
ditions can also be satisfied by properly choosing the free 
parameters. 

Finally, we note that da Silva argued, in the HMT 
setup, that the coupling constant A can still be different 
from one [12[. If this is indeed the case, then one might 
be concerned with the strong coupling problem found in 
other versions of the HL theory @, [23, [28|. However, 
since the spin-0 graviton is eliminated completely here, 
as shown explicitly in [To|, [ID, UU, this question is au- 
tomatically solved in the HMT setup even with A ^ 1 
[29l | . It should be noted that da Silva considered only 
perturbations of the case with detailed balance condi- 
tion, and found that the spin-0 mode is not propagating. 
It is not clear if it is also true for the case without de- 
tailed balance. The problem certainly deserves further 
investigations. 

Note Added: A preprint 25!] appeared in arXiv al- 
most simultaneously with ours. These authors also stud- 
ied spherically symmetric static vacuum solutions similar 
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to those in Sec. IV, but did not consider the subjects pre- 
sented in the other parts of this paper. 
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+ ^(^ + vv)-^(^ + vV)}, 

F H,3) = \ {2V fc V^ )fc - V 2 /;-? - (V fc V,/*') g ij } , 

(A.2) 

where 



Appendix A: Functions (-Fs)y and F^ » 



p{ [2K l: > + V J VV) - ^[2K + V 2 (p)g* 3 



(A.3) 



The geometric 3-tensors F u and n ~. defined in 



Eq. (|2.14j ) are given by 
1 



%,2) 



2 9ij, 

1 

2 (J2y - ViVj) R - l -g l3 (R - 4V 2 ) R, 



(F 3 ) i? = V 2 i? l3 - (V,Vj - 3i^) R - 4 (i? 2 V 



"2^3 (3RkiR kl + V 2 R-2R 2 ) 



3 (iiy - ViVj) R 2 - -gij (R - 6V 2 ) i? 2 , 
(iZy + ViV,-) (R k iR kl ) + 2R (R 2 ).. 
+V 2 {RR l3 ) - V fe [V, (ilflj-k) + Vj (ARifc)] 
-i 5i ,[(i?-2V 2 ) (irf) 



(^) l7 = 3(i? 



2V fe Vi (flfl fci )] 
3 
2 



V 2 fi? 2 ) 



V fe ( Vi (i? 2 ) 



1 



R k R l m R™ - 3V fc V; (i? 2 ) 

72 



A'/ 



2ViVj (V 2 i?) - 2 (V 2 i?) i?„ 
1 



(ViR) (VjR) 



2' 

V 4 R tJ - V fc (V.V 2 ^ + V 3 - V 2 R k ) 
- (V<itf) (VjRi) - 2 (V fe ^) (Vfc^,) 



(Vi?) 2 + 4V 4 i? 



"2 5ij 



(A.l) 



(V fc i?, m ) 2 - 2 (V fc V;V 2 # H ) 
- <p { (W + V 2 <p) R l > - 2 + V J V fc </>) R 
- 2(2x1 + V 4 V fc ^)i? jfc 



- ( 2A 9 - R ) { 2K tJ 



Appendix B: Function FV n -. in Spherical Static 
Spacetimes 

In the spherically symmetric static spacetimes de- 
scribed by the metric (|3.1[) . the function , defined 
by Eq. (|A.2[) is given by 



[ F ( v ,2)] 



-2v 



! "(l-A 9 r 2 ) - (1 + tV) 
(A g r 2 - l) + (3 + rv 1 ) 



i/ip' 



e 2 »{A g r 2 -l) +1 



- 2e^ +v 

+ 2e» +v {2-r t x')v'Y r i 8 T j 
rv' ~{\-e 2v 



1 



o / 2v \ I I oA 2ft 

6 — rv — e )v<p — 2A g re tp 



J „2v \ ,,/ 



- 2e^ +v {2 - rv -e"J/i' 
+ 4e" + "(i/- A g re 2 ») W,, 



e 2 ^V 2 -i) + i 



ipip 



ij 2r 2 

+ (cp' + tpv'^j iff + 2e» +v (V + <pA 
+ 4spef i+v (t/ - A 



2A s r^e 2 V 



+^ e - 4v ^(u'-A g re 2u ) ( p'' 

+ rtp(ip' + 2e» + ^v" 
r(i/- A g re 2u ^<p' 2 
hip' + 4e^ +,y ) v 



+ ri 
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F 



(V,3)J 



-(w - 2re^ +u - A g r 2 ipe 2u ^v'ip' 
-2e^ v (A g r 2 e 2v tp' + ipv'^j 

-2v r 

— ripv'ip" + <p' 2 - ip{rv' + 2)u'(p' 



8 r , 5 r , 

L J 



+ 2e^ +u U + ripv'/j,' - 2</jz/) 

-np(y + 2e" +l, y" 

+ &r[tpv' — ip J + tp v ip 

+ e» +v Ur<p(v' - n')v' + 2ipv' 

+ 2r(fi' -3v')ip' (B.l) 

Appendix C: Impulsive Parts of , 

From Eas.(|7.ip and (|7.4p we find that [-F^.n)]^ given 
by Eq. (|B.l[) takes the form, 



F 



S(r-r ), (C.l) 



where 



i 7 , 



(V.3) 



-2i/ 



1 - A r 2 ) - (1 + ri/) 



"2^ 



r;/ — 1 — e 



-2y 



2r 2 



e 2 ^A 3 r 2 -l) 4 I 



-4// 



(V-A ff re 2 ") fo/] 



+ (V + 2e^+^ [*/] 



-—re 
2 



2(^-^' + e^ +l/ ) fo/]" 



£V (c.2) 
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